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Q-r ABSTRACT 



We calculate spectral energy distributions (SEDs) of steady accretion discs at high 
accretion rates, as appropriate for bright QSOs, under the assumption that the outer 
parts are heated sufficiently to maintain marginal gravitational stability, presumably 
by massive stars formed within the disc. The SED is independent of the nature of these 
auxiliary sources if their inputs arc completely thermalizcd. Standard assumptions are 

■ made for angular momentum transport, with an alpha parameter less than unity. With 
, these prescriptions, the luminosity of the disc is sensitive to its opacity, in contrast to 

standard discs powered by release of orbital energy alone. Compared to the latter, our 
discs have a broader SED, with a second peak in the near-infrared that is energetically 
comparable to the blue bump. The energy in the second peak increases with the 
£NJ \ outer radius of the disc, provided that the accretion rate is constant with radius. By 

■ comparing our computed SEDs with observed ones, we limit the outer radius of the 
^ (— | disc to be less than 10 5 Schwarzschild radii (Rs), or about one parsec, in a typical 

fS ^ QSO. We also discuss some properties of our minimum-Q discs in the regions where 

auxiliary heating is dominant (10 3 — 10 5 i?s). 
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1 INTRODUCTION 

The standard theoretical paradigm for the central engine in quasars and their radio-quiet kin (QSOs) is a viscous accretion 



■ disc surrounding a massive black hole. Direct evidence for accretion discs, such as double-peaked emission lines (Eracleous 



fe Halpern L994), is important to seek but hard to come by. Perhaps the best reasons for belief in this paradigm are basic 
considerations of energy and angular momentum. Such discs are the most plausible astrophysical mechanism for converting 
rest mass to radiation with high efficiency, and if the relics of QSOs reside in galactic nuclei, then efficiencies > 10% are 



required (soltan 1982; Chokshi & Turner 1992; Yu & Tremainc 2002). On the other hand, the maximum specific angular 



momentum of a Kerr black hole, GM/c « 1.4Afgkm s _1 , is far less than that of most mass in galaxies. (Ms = M/W 8 Mq, 
where M is the mass of the black hole.) Accreting gas must be separated from its angular momentum, and a viscous disc is 
the most natural mechanism. 

Most of the binding energy of the disc is in its inner parts, so that the outer radius of a standard accretion disc is almost 
irrelevant to its bolometric luminosity. In contrast, most of the angular momentum is in the outer regions, so that the outer 



radius of the disc is sensitive to the initial angular momentum of the gas that feeds it. In a previous paper (Goodman 2002, 
henceforth Paper I), it was argued that luminosity and angular momentum are coupled by the requirement that QSO discs 
be stable against their own self-gravity. 



It is well known that gravitational stability is problematic in the outer parts of QSO discs (Shlosman & Begelman 1987). 
The threat to the standard paradigm is that a strongly selfgravitating disc is likely to fragment completely into stars, leaving 
insufficient gas to fuel the QSO. As discussed in Paper I, viable solutions to this difficulty fall into several categories: 

(i) enhanced angular momentum transport, not necessarily by viscous processes but at rates corresponding to a viscosity 
parameter a>l; 

(ii) auxiliary heating in excess of what is provided by dissipating orbital energy, so as to reduce the gas density and 
selfgravity of the disc; 

(iii) replacement of the outer disc with a very dense star cluster, whose collisional debris supply a disc of small radius and 
negligible selfgravity; 
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(iv) relatively low initial angular momentum for the gas, which therefore circularizes at small enough radius so as to avoid 
selfgravity. 

A number of options in the first category were briefly considered in Paper I, including accretion driven by bars or global spiral 



waves (e.g. Shlosman & Begelman 1989) or by magnetized winds (e.g. Blandford & Payne 1982). But it was argued that none 



of these options is likely to achieve a supersonic accretion speed, and hence that discs cannot be stabilized much beyond one 
parsec by any of these mechanisms alone. The third category was tentatively rejected on the grounds that remnants of such 
star clusters are not observed in present day galactic nuclei. Recently, Pariev et al. ( 2002 ) have proposed discs whose thickness 



is supported primarily by magnetic pressure. If this is possible, it would imply lower disc densities on average, but the field 



might well squeeze the gas into dense clumps and thereby actually exacerbate selfgravity. At any rate, Pariev et al. (2002) do 



not apply their model to the cool outer regions beyond 10 -Rs- 

The present paper will focus on some implications of the second category of solutions, which appear to be a natural 
compromise between a purely gaseous disc and a star cluster. It is reasonable to suppose that part of the gas fragments into 
stars, and that the energy released by nuclear fusion and other stellar proc esses (supernovae, ste llar-mass black holes) may 
sufficiently heat the rest of the gas so as to prevent complete fragmentation ( Collin fc Zahn 1999a , b|) . There is a great deal of 
evidence that such a feedback cycle operates in the discs of spiral galaxies on kiloparsec scales. 

Observed SEDs of typical quasars differ markedly from classical theoretical predi ctions in whi ch the disc is assumed to be 
geometrically thin, optically thick, steady, and heated solely by viscous dissipation (Pringle 1981). To a first approximation, 
the typical SED is flat in a \F\ plot over many decades in wavelength (Elvis et al. 1994). Relative to the classical predictions, 
there is excess emission at both X-ray and infrared wavelengths. The former is conventionally ascribed to comptonization in 
a hot corona at small radii ( Shapiro et al. 1976| ) , and the latter to passive reprocessing in warped or flared outer parts of the 
disc (Sanders et al. 198£). 

We suppose that the infrared excess may be due to the energy inputs required to stabilize the outer disc against its own 
selfgravity. We obtain a lower bound on the auxiliary heating needed to stabilize the disc for given values of the macroscopic 
parameters: namely, the black hole mass (M), accretion rate (M), and disc outer radius (r max ). As shown in Paper I, and 
confirmed here with more realistic opacities, the inputs required for gravitational stability increase with the outer radius of 
the disc, r max . Paper I argued for an upper limit to r max based on the energy available from plausible sources, such as fusion 
or accretion onto stellar-mass black holes. In this paper, we find limits to r max from the SED. 

It is not obvious that the auxiliary inputs should be completely thermalized, but if we assume this, then the SED can 
be predicted. The observed infrared emission of typical QSOs may be due largely to reprocessing of light emitted from the 
inner parts of the disc, as conventionally supposed. But by attributing all of the infrared light to the auxiliary energy inputs, 
and insisting that these inputs be sufficient for gravitational stability, we obtain bounds for r max . These bounds depend upon 
other parameters, especially the mass of the black hole, the accretion rate, and the viscosity parameter a. We explore these 
dependencies. The meaning of r max constrained by this method is the radius within which M is sensibly constant; obviously 
the disc can be extended indefinitely if the accretion rate and mass at large radii are sufficiently small. Hence our limits on 
7" max are best translated into upper limits on the initial angular momentum of the gas that is accreted. 

The outline of our paper is as follows. §2 lays out the physical assumptions and governing equations for our disc models. 
These are the same as for the classical steady thin disc, with the one important exception that wherever the classical model 
would be gravitationally unstable, we invoke just enough auxiliary heating to stabilize it. A discussion of opacities becomes 
critically important, because the requirement of marginal stability fixes the density and temperature at the midplane (for 
given M and a); the flux escaping from the disc, and hence the amount of auxiliary heating needed, then depend upon the 
optical depth. The computed radial structure of the disc is presented in §3 for parameters representative of bright quasars. 
Since our assumptions are no different from the classical ones at small radii, we emphasize the properties of our discs in the 
marginally selfgravitating region r > 10 3 R S ~ 10" 2 pc, and we compare our SEDs with those presented by pivis et ah] (|1994| ) 
to obtain upper limits on the outer radius of the disc and the initial angular momentum of the gas. In the final section, we 
summarize our conclusions, issue the necessary caveats, and discuss directions for future research. 



2 THE SELFGRAVITATING ACCRETION DISC MODEL 

We start by considering the standard steady a disc. In this model, gas is accreted onto a central black hole of mass M ~ Mg 
in a steady, Keplerian, geometrically thin accretion disc at a constant rate M. Accretion is driven by viscous mechanisms, 
but the exact mechanism (probably magnetorotationally-driven turbulence) is unimportant here. Newtonian equations will be 
used. These are not accurate near the inner edge of the disc, but the effects of selfgravity are at large radii where relativistic 
corrections are unimportant. 

We assume that each annulus of the disc radiates as a blackbody with temperature T e g(r). When accretion is the only 
source of thermal energy, 
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aT 4 ff = ^-M'n 2 , (1) 
on 

where Q = (G"A//r 3 )~ 1/2 , and we have defined for notational simplicity M' = Mil — \/r sa i n /r) where r m i n is the inner radius 
of the accretion disc. A zero-torque boundary condition has been assumed. In a relativistic thin-disc analysis, r m i n would be 
the radius of the marginally stable orbit. In our newtonian framework, we set r m i n = i?s/4e, where Rs = 2GM/c 2 is the 
Schwarzschild radius and e ~ 0.1 is the assumed radiative efficiency of accretion. 

An a prescription for the viscosity is adopted, v = ac s hf3 b , where c s = \ppfp is the isothermal sound speed at the disc 
midplane, h — c s /Q. the half thickness, and (3 = J5 gas /p is the ratio of gas pressure to total pressure at the midplane. The 
viscosity is related to the accretion rate by 

where E = f^°° pdz is the surface density at a given radius. The parameter b is a switch that can be either 1, so that 
viscosity is proportional to gas pressure, or 0, so that viscosity is proportional to total pressure. Magnetic pressure is assumed 
to contribute negligibly to the thickness of the disc, even if magnetic stresses dominate angular momentum transport. 

Because it will turn out that the optically thick assumption fails at large radii, we now turn to a study of the relationship 
between the midplane temperature T and T a a. We assume radiative energy transport in the vertical (z) direction and write 
r — Kpdz for the optical depth at the midplane. In the optically thick limit t>1, the diffusive approximation is applicable 
so that T 4 ~ i~T 4 B , with a numerical factor of order unity that depends upon the vertical variation of opacity « and the vertical 
distribution of viscous dissipation. The latter is not yet predictable from theory, despite major advances in understanding 
magnetohydrodynamic disc turbulence. If one assumes constant dissipation per unit optical depth, then 

where the "+1/2" follows from the Eddington approximation to the boundary condition at the photosphere. In the opposite 
limit of low optical depth, the radiative flux from one side of the disc scales as crT 4 ff ~ raT 4 . Again a numerical factor occurs. 
One should distinguish between absorption and scattering opacity, but at the large accretion rates of interest to us, optically 
thin conditions tend to occur where the opacity is primarily absorptive. For a vertically isothermal disc, aT 4 s = 4t(tT 4 when 
r < 1. Simple interpolation between the optically thin and optically thick limits yields eq. (^) below, which satisfies the 
physical requirement that T e g ^ T for any r (in fact, the minimum of T 4 /T 4 S is 1.112 and occurs when r = v/2/3). 

We represent the effective radiation pressure by eq. ([?]). In the optically thick case, in view of eq. (^|), this reduces to the 
usual p rac i = aT 4 /3 relation, but in the optically thin case, p ra d — ► 2r 2 aT /c: one power of r represents the inefficiency with 
which photons are radiated, and the second is the fraction of these photons that transfer their momentum to the gas, thereby 
helping to support the thickness of the disc. The numerical factor is correct for an isothermal slab of midplane optical depth 

T < 1. 

We also must know the opacity in the disc (equation ^). We use the opacity tables of Iglcsias & Rogers (1996) for high 



temperatures and those of Alexander & Ferguson (1994) for low temperatures, with X — .70, Z = .03. Figure |fj illustrates 



the opacity as a function of density and temperature. The two opacity tables overlap in the region 3.75 < logT < 4.1, and 
Figure ^ demonstrates their agreement. Additionally, the loci of two disc models are projected onto the plot. Extrapolation of 
the opacity tables in the low density regime is necessary. This is done by simply assigning k its value at the low-density edge 
of the opacity table (using the same value of T). This extrapolation is presumably correct for high temperatures (T > 10 4 K) 
where the opacity is dominated by electron scattering. At lower temperatures, the extrapolation is less certain, but we feel it 
is reasonable as a first approximation. Also, for logT < 3, ft is assigned its value at the low-temperature edge of the opacity 
table, which is approximately 10 0,76 cm 2 g _1 . 

In summary, our governing equations for the disc are as follows (we take m — .62 tuh as the mean molecular mass): 

^Tcff = ^-M'n 2 (3) 



,8 2 At 



(5) 



P c s E = (6) 

i-ROl 

Prad = -^ T cS ( 7 ) 

pkT ^ 
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as a function of density and temperature, as provided by Iglesias & Rogers (1996) (high-T region) 
i) (low-T region) for X = .70, Z = .03. Contours represent integer values of the logarithm of the 



dashed (6 
r m ax = 10 5 R S 



). The overlap region 3.75 < logT < 4.1 is plotted with contours from both sources; their agreement is illustrated. The 
0) and dotted (b = 1) lines show trajectories of (p, T) for typical examples of our disc models (Figure ^ for b = 0), up to 



p = — Be — (9) 

+ Prad 

E = 2ph (10) 

c 2 = Pgas+p r ad (12) 
P 

K = K(p,T) (13) 

For an accretion disc with given parameters M, M, a, and b, we can solve these eleven equations for the eleven unknowns: 
T c ff, T, r, E, /3, c s , Prad, Pgas , p, h, and k at every radius r, which is related simply to Q. 

The outer parts of the accretion disc will be prone to selfgravity if Toomre's stability parameter, 

is < 1. In our model, we assume there is some feedback mechanism, perhaps star formation, that supplies just enough additional 
heat to the disc as to prevent Q falling below a minimum value Q m i n ~ 1. The energy source for these auxiliary inputs must 
be something other than the orbital energy of the gas, for example, nuclear fusion in stars. The quantity of these inputs is 
fixed by the properties of the disc: in particular the opacities, since the more easily radiation escapes, the more the disc must 
be heated to maintain stability. The nature of the auxiliary sources is not important to the spectral energy distribution as 
long as their inputs are completely thermalized. 
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When Q = Q m i n , eq. (g), which assumes the disc to be heated by release of orbital energy only, is no longer appropriate; 
we replace it with 

p= ^Fn- (15) 

and T e ff must be obtained from equation ^. 

Once T c fi is known as a function of r, we compute a theoretical spectral energy distribution for the disc as a weighted 
sum of Planck functions: 

LX = ^ r L„ 27Vrdr eMhc/XkT cB (r)]-l (16) 

This formula does not allow for differences between effective and col or temperature. Suc h differences are most likely in those 
parts of the disc whose emission is dominated by passive reprocessing ( [ganders et al. 1989) or are optically thin (to absorption). 
In fact we find that much of the auxiliary heat input occurs in regions that are marginally optically thick, since these regions 
cool most efficiently. We feel that eq. (|l^) is adequate for a first investigation: it may somewhat overestimate the wavelength 
of the infrared peak associated with the auxiliary heating, but the area under the peak should be roughly correct. 



3 RESULTS 

We define Ze = Lo/Le, where Lq = eMc 2 is the luminosity corresponding to a nonselfgravitating (i.e., with no auxiliary 
heating) disc with accretion rate M , and Le is the Eddington luminosity. In our models we always set e = .1 and Q m in = 1- 

The radial structure of a typical accretion disc (Ms = 1, Ze = -5, a = .01) is presented in Figure Q for 6 = 0. The outer 
radius of r max = 10 7 i?s will be shown to be unrealistic, but this model illustrates several key regimes in the disc. 

The inner region (2.5 <r/Rs< 10 2 ) is radiation dominated (/3 -C 1) and the temperature is high enough to ensure that 
the opacity is dominated by electron scattering, hence constant. In view of equation |l], T e g oc r~ 3 ^ 4 . When Q would otherwise 
drop below 1 (here, at r « 10 3 i?s), we hold Q constant, which ensures that p oc r -3 . In this auxiliary heating region, if 
K = constant, then T e g oc r _3//8 . After the disc becomes optically thin (r < 1) we find multiple solutions. The one with the 
least optical depth corresponds to the continuous line in all panels; it reaches constant T and k. There is an optically thick 
solution with T oc r -3 / 4 . Since the region in which we find multiple solutions is so dependent on the extrapolation of opacity 
into low-density and low-temperature regimes, the interpretation must be approached with caution. However, the presence 



of multiple solutions raises the possibility of thermal instability and time-dependence (see, for example, siemiginowska et al 



199C ) . Note that in our solutions, M is constrained to be constant throughout the disc. Several of the panels in Fig. |2j contain 
only one curve in the region of multiple solutions (r > 10 pc): this is because p and c s are uniquely determined by r, M, a, 
and Q via eqs. ( [t5| ) and ^ when 6 = 0. 

For 6=1, the effective temperature has the same T e g oc r" 3 / 4 behavior in the inner part of the disc, but differs from the 
6 = case in the auxiliary heating region. For n oc constant in this region, in the optically thick solution, T oc r -1 ' a & T e g oc 
r -1 ' 4 , and in the optically thin case, T oc const. & T e g oc r~ 3 ^ s . 

The calculated SED of this typical model is presented in Figure ^, for outer radii of r max = 10 3 , 10 4 , 10 s , 10 6 , and 
10 7 i? s . (All of our SEDs assume that the disc is viewed face-on.) Since T e g may have multiple solutions in certain regions 
of the disc, the SEDs are plotted as bounded regions with the top (bottom) boundary corresponding to the largest (least) 



T e ff. Also plotted is the arbitrarily normalized mean energy distribution (MED) from 47 observed quasars (Elvis et al. 1994). 
The authors caution that the dispersion about this mean is large, of order "one decade for both the infrared and ultraviolet 
components when the MED is normalized at [A = 1.25 fxm]." The r max = 10 3 7?s case has an expected A/a oc A -4 / 3 behavior 



since selfgravity is not an issue. Relativistic effects (Laor & Netzer 1989; Sun & Malkan 1989) only influence the SED at and 



blueward of the short wavelength peak, since this emission conies from the innermost parts of the disc. 

Figure |j| shows the dependence of the infrared bump on the mass of the central black hole. (Note that the wavelength 
corresponding to the blue bump has an expected M 1 ' 4 dependence.) The Eddington ratio as defined above is kept constant 
at Ze = -5: note that this parameter defines M rather than the total luminosity of the disc, which includes the auxiliary heat 
inputs. In the top two panels, r max is set to 10 5 i?s, and the relative prominence of the infrared bump grows with M because 
the reduced radius at which Q m i n = 1 decreases with M (Paper I). The bottom two panels have r max fixed at 1 pc; here larger 
values of M tend to reduce the infrared bump. This dependence is illustrated in Figure [?]. Solutions for the disc are considered 
incompatib le wi th the observations if A/a at the maximum of the infrared bump is more that ten times that of the blue bump. 

Figure |5(a)| shows that while the blue bump increases with mass accretion rate, the energy output of the infrared bump 
is hardly affected, for the 6 = case. This invariance is caused by two competing factors as M is increased: T e ff is higher 
in the selfgravitating region, but because T is also higher, the onset of the optically thin regime, and thus a sharp incre ase 
in T c ff as shown in Figure [| is delayed. Similar considerations explain the behavior of the infrared bump in Figure 5(b) for 
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Figure 2. Dependence of an assortment of parameters on radius for our b = (i.e., viscosity oc total pressure) canonical case of Mg = 1, 
£e = -5, a = .01, and Q m j n = 1. All units are in cgs except for those represented as reduced variables. The integrated disc mass 
J r r 27rr'Sdr' is in units of the central black hole mass M, but the most unstable mass c 4 /G 2 E is in units of solar masses. The SED 
for this model is the vertically-hatched region in Figure ?(a). 



6=1. For low accretion rates, the midplane temperature drops low enough that the disc becomes optically thin before the 
end of the disc (at r max = 10 5 i?s), which drives up T e ff . 

In Figure |(| the same argument explains the behavior of the infrared bump in both the 6 = and 6=1 cases. The onset 
of the optically thin selfgravitating region is at smaller radii with increasing a for both cases, but T c g is higher in the optically 
thick selfgravitating region for 6 = 1. 

It is clear that additional heating in the outer parts of our accretion discs results in a second bump in the SED. Assuming 
that this outer-disc contribution to the SED can be no more prominent than infrared emission in the observed MED, we can 
already place tight constraints on the outer radius of these disc models: r max < 10 4 — 10 5 i?s ~ 0.1 — 1 pc. 



4 SUMMARY AND DISCUSSION 

We have estimated spectral energy distributions (SEDs) of bright QSOs using standard assumptions, with one addition: 
where the disc would otherwise be gravitationally unstable, we have postulated additional sources of heat, other than release 
of orbital energy by accretion, just sufficient to maintain gravitational stability. These sources become necessary beyond 
~ 10' ! i?s for typical parameters. Assuming their energy inputs to be locally and completely thermalized, we have calculated 
their contribution to the SED and luminosity of the disc, which occurs primarily in the red and near infrared. The larger the 
disc, the more auxiliary heating is required. For typical black-hole masses and accretion rates inferred from the blue bump, 
the auxiliary inputs actually exceed the power derived from accretion if the disc extends beyond 10 4 — 10 5 i?s, or about one 
parsec. This would be incompatible with the typical SED of bright QSOs, which is approximately fiat in XF\. 

Paper I placed similar limits on r max from energetic arguments. It was assumed that the stars or small black holes that 
heat the disc also form within the disc, and that the mass in these objects is at most comparable to that of the disc. The 
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Figure 3. Calculated spectral energy distributions for our accretion disc models with Mg = 1, £g = .5, a. = .01. The five SEDs, 
in order of increasing total luminosity, correspond to the five values of r max /R s = 10 3 , 10 4 , 10 5 , 10 6 (diagonally-hatched), and 10 7 
(vertically-hatched). Multiple solutions for T e g (r) are responsible for the ambiguity of the SEDs plotted as bounded regions. 



present limits do not rely on these assumptions. That is to say, if the disc were heated by a much larger mass in stars, the 
limits of Paper I could be relaxed, but those based on the SED would still apply. 

It has been assumed that the accretion rate is constant with radius, so that M at large radii can be derived from the 
luminosity in the blue bump (which has no contribution from the auxiliary sources). If M at r > 1 pc is several orders of 
magnitude less than it is at r < lO^-Rs, then the disc could be much more extensive than we have supposed. The Stefan- 
Boltzmann law implies that the observed luminosity at 10 fj, must in any case come from r > 2(ALa/10 46 erg s _1 ) 1//2 pc, but 
this does not require that the surface density at that distance is as high as in a constant-M disc. Therefore, it is useful to 
rephrase the limits in terms of the initial angular momentum of the gas supplied to the disc. At least on a time average, M 
should be constant inside the radius corresponding to the initial angular momentum. The relationship is 

^ 1/2 pekm s _1 . 



Jo ~ 660 (r pc M 8 ) 



(17) 



This is quite small compared with the product of virial velocity (~ 300 km s _1 ) and scale size (~ lkpc) of QSO hosts, and it 
may be important to ask where gas with such low angular momentum comes from. Furthermore, the mass of a gravitationally 
stable disc that obeys our constraints on the SED is generally much less than that of the black hole, so that the disc must be 
replenished many times over to grow the black hole by accretion. 

Our treatment of the disc is highly simplified and certainly crude compared with many past efforts. We have ignored 
relativistic effects, adopted a one-zone model for the vertical structure, used Rosseland mean opacities without distinguishing 
between scattering and absorption, taken a constant molecular weight, and represented angular- momentum transport by the 
usual viscous a prescription. We feel that these simplifications are justified by the strong dependence of the auxiliary inputs 
and the SED on the outer radius of the disc. A more detailed treatment of the physics that one actually understands seems 
unlikely to change our conclusions concerning r max - Radical enhancements in transport (equivalent to a > 1), or magnetic 
pressures S> p ga a (as suggested by Pariev et al, 2002) could make some difference, but these are not yet understood. 



Despite what has just been said, a more elaborate treatment of vertical structure and radiative transfer might point to a 
redistribution of the auxiliary energy inputs in wavelength, if not their total contribution to the disc luminosity. Unfortunately, 
one has no reliable predictions for the vertical distribution of purely viscous heating, much less of the auxiliary sources 
postulated here, which will limit the credibility of detailed vertical models. 

Even if the actively accreting parts of QSO discs are smaller than 0.1 pc, it is still possible that selfgravity is important 
in them, and therefore that they form stars. Quiescent galactic nuclei with black holes, even in early-type galaxies, often show 

120011; 
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kinematic evidence of compact stellar discs ( Glebhardt et al 
(albeit on scales ~ 10 2 pc) are usually accompanied by AGN activity ( Sandere 



199£ 



le Zeeuw et al,_ 

1999^ 



Heckman 



2002). Nuclear starbursts 



The black hole in our 



own Galaxy, though not an AGN and estimated to ha ve a very low accre tion rate ( Quataert et al. 1990 , is surrounded by 
what appear to be young high- mass stars at r < 0.1 pc ( Krabbe et al. 1995 ). More theoretical attention should be paid to star 
formation in these extreme environments, where the densities, temperatures, and tidal fields are much higher than in normal 
giant molecular clouds. 
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(c) 6 = 0, r max = 1 pc (d) b = 1, r max = 1 pc 

Figure 4. SEDs for the models with = -5, a = .01. The four SEDs shown here, in order of increasing blue bump, correspond to 
logAfg = —2, —1, 0, 1. The outer radius is set to r^ax = 10 5 Rg in the top two panels, but it is fixed at r max = lpc in the bottom two 
panels. Compare the bottom two panels to Figure M. 

We thank Iskra Strateva and Jonathan Tan for helpful discussions. 
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Figure 5. SEDs for the models with Mg = 1, a = .01, r max = 10 5 -Rg. The five SEDs shown here, in order of increasing blue bump, 
correspond to = .01, .03, .1, .5, 1. 
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Figure 6. SEDs for the models with Mg = 1, Zg = .5, r max = 10 5 i?g. The three SEDs shown here, in order of decreasing (increasing) 
infrared bump for b = (b = 1), correspond to a = .01, .1, .3. 
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Figure 7. Locus of allowable (M. r ma x), with Ze = -5 and a = .01 fixed. The bottom dashed (dotted) contour represents solutions for 6 = 
(6 = 1) for which (A/ a )ir bum P /(A/A)bluc bump — I? and the top contour represents solutions for which (A/a)ir bump/ (A/a )bluc bump — 
10. Thus, the diagonally-hatched region is the obscrvationally-constrained approximate forbidden region for 6 = 0, and the vertically- 
hatched region is the forbidden region for 6 = 1. In both cases r max is forbidden to be much larger than ~ lpc. 
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